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ABSTRACT 



The basic Mitrovic*s method is an effective technique in the 
analysis and design of linear feedback control systems. Mitrovic’s 
method has been successfully applied to analyze feedback control systems 
with single nonlinearities. The objective of this work was to employ 
Mitrovic’s method, which permits the variation of two coefficients of 
a characteristic equation, in the analysis of a control system with two 
gain-variable nonlinear feedback paths. 

After predictions of system performance were made, the predictions 
were tested by simulating the feedback control system on a Donner 
Scientific Corporation analog computer. Model 3100. Computer results 
were analysed in order to support or reject prediction techniques. 
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CHAPTfclR I 



USE OF THE BASIC MITROVIC'S METHOD 
In general, the analysis and design of feedback control systems 
centers around the solution to an ordinary, linear, differential equation 
with constant coefficients. Specifically, the application of LaPlace 
transforms to the differential equation produces a transfer function 
which is a ratio of the output signal to the input signal. By proper 
placement of the poles and zeros of the transfer function, the frequency 
and/or time response are adjusted to produce results which meet with some 
set of specifications for the system. This procedure may be followed 
because the form of the solution to a linear differential equation is 
unique. The resulting solution is predictable and invariant. 

The introduction of a single nonlinear element in a control system 
produces a nonlinear differential equation where one or more of the co- 
efficients may be variable. There are as many specific solutions to a 
nonlinear differential equation with variable coefficients as there are 
values of the coefficients. Therefore, the accurate prediction of system 
performance is stringently curtailed. The use of describing functions 
provides an adequate solution to the problem of one nonlinear element 
within the limitations imposed on the describing function itself. 

Consider a feedback control system with a transfer function GCjCJ ) 
and a nonlinear device which may be represented by a describing function 
( j 6) ) both in the forward path as sho\>m in Fig. I-la. Then the 
transfer function for this system is: 

1 Gu(3a>') 



1 



• 

1 

^ . 



I 




' ♦ 




f 

i ^ 
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from which the characteristic equation may be found by setting the 
denominator of Eq. (I~l) equal to zero. Thus, 



1 G GjjCjw') = 0. 



( 1 - 2 ) 



The solution to Eq. (1-2) is: 






Gob") 



(1-3) 



Where G^(j6l) ) is a gain nonlinearity only and has no phase angle 
associated \;ith it - which is the case for such nonlinearities as satura- 
tion, ideal relay, dead zone, and relay with dead zone - the problem is 
conveniently solved on a gain-phase plot. This graphical solution is 
illustrated in Fig. I-lb. The system operates in a limit cycle at the 
intersection of the two curves. This is shown as point "A” in the 
figure. A root-locus plot may be used to arrive at the same result 
providing the single nonlinear gain is not a function of frequency. 

When two nonlinear devices appear in a control system, only a 
family of gain-phase plots or root-loci plots can provide a complete 
solution to the problem. Thus the problems of analysis and design become 
complex and cumbersome. Analysis would be convenient if: 

1. Some method of simultaneously looking at i\7o nonlinear 
gains were available, and, 

2. Some relationship between these two nonlinear gains could 
be found. 

Mitrovic's method [1] provides the answer to the first requirement for 
control systems in which each nonlinear gain appears in a different co- 
efficient of the characteristic equation. The purpose of this thesis 
is to fulfill the second requirement for a particular feedback control 



3 



system. 



The mathematical derivation of Mitrovic*s method is based on a 
theorum by Cauchy and will not be pursued here. The algebraic manipula- 
tions give some insight into Mitrovic*s method and will serve to provide 
the reader with some feeling for the method. 

Consider a feedback control system with a transfer function 



G(s) 



K 

CS--PJ 



The characteristic equation for this system is 



(1-4) 



* pip2,s * K = 0 , 

which is of the form 



(1-5) 



Aj Ai S = 0 . (1-6) 

If the coefficients to be varied are and then the characteris- 

tic equation may be rewritten as 

AjS’ -► AjS* t B.S - = 0 , a-7) 

where the values of s have the general form of 

5 = -V a-s) 

and lie in the second quadrant of the s-plane. Then the values of s 
which must be substituted in Eq. (1-7) are: 

-1) - j 

Substituting the values of s from Eqs. (1-8) and (1-9) in Eq. (1-7) 

and requiring the real and imaginary parts to go to zero independently 

4 



/ 




yields : 



* Bo = 0 (I-IO) 

io„f^[-A3«^(l-4f) - k^Z\u^ * BJ =0. »-”) 

The solution of these two simultaneous equations, after first dividing 
Eq. (I-ll) by becomes : 

Bo = -[AsW* it)) - (1) A;,60^] a-12) 

Bi = [Ajto‘(l-4's’') ♦ AiWjz's')]. 

By plotting the results of Eq. (1-12) on a coordinate system for 

constant ^ lines as varied, it is possible to map radial lines 

from the s-plane onto the plane. Thus the result of simultan- 

eously varying two coefficients of the characteristic equation is plainly 
evident from a single plot. l\Jhen the ^ ^ line (s = in the s-plane) 

is selected as radial line to be mapped, the stability curve is plotted 
on the Bq - plane. 

Mitrovic's method may be expanded to solve any order characteristic 
equation for any two coefficients. The generalized solutions for Mitrovic*s 
equation pairs are listed in Appendix B, Table I. In a generalized solu- 
tion, certain functions of only ^ repeated, such as 2^ 

in Eqs. (1-12). A listing of these functions of ^ is made in Appendix 
B, Table II. 

The analysis of feedback control systems with single gain variable 
nonlinearities using Mitrovic*s method is straight forward and produces 
results consistent with root-locus, describing function, and analog 
computer techniques. Since the author had no previous experience in 



5 



either Mitrovic*s method or analysis of feedback control systems with 
two gain variable nonlinearities, throe control systems were selected 
and analyzed before study on the two gain-variable nonlinear problem 
was begun. The results of the study of these three systems appear in 
Appendix A and corroborate the results obtained by LT P. L. WILSON, USN. 
[ 6 ]. 
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CHAPTER II 



THE FEEDBACK CONTROL SYSTEM WITH ACCELERATION AND VELOCITY 

FEEDBACK CCMPENSATION 

A block diagram of the system to be studied in this thesis is 
shown in Fig. II-l. The root locus of the transfer function for the 
system, without acceleration or velocity feedback paths, is plotted on 
the s-plane in Fig. II-2a. For a gain of 60, the system has complex 
roots in the right half plane and is unstable. The addition of velocity 
and acceleration feedback compensation, when saturation is not present, 
produces a characteristic equation which is 

+ 63s^ + 62s + 60 = 0, (II-l) 



when K and K are both unity, 
at 

Eq. (Il-l) may be factored to give 



( s + 62.016) ( s + 0.984s + 0.984) = 0. 



(n-2) 



Thus, the system has been compensated to approximate a second order 
system with a damping ratio, ^ , equal to 0.5 and a natural frequency, 
Co ^ > equal to 0.984. The root locations for the compensated system 
are shown in Fig. II -2b. In response to a step input, the system would 
have a maximum overshoot, 1.15 and a settling time of about 

four seconds. 

The addition of two nonlinear elements produces a characteristic 
equation with two nonlinear coefficients. This equation is: 



+ (3 + 60N,K )s^ + (2 + 60N K )s + 60 = 0 , 
^ a 1 L 



(n-3) 
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FIGURE II- 1 

■r 

Feedback control system block diagram. 



W ■■ 




where and represent the instantaneous gains of the velocity 

and acceleration saturating amplifiers, respectively. The instantaneous 
value of Nj and N 2 may be either unity or a magnitude determined by 
Eqs. (II-4). 



N = 
1 



"2 = 



sat 



E 



sat 



e ) 



(II-4) 



The root locations for Eq. (II-3) vary with time as 0 and 0 vary. 

This is evident if the block diagram is first reduced by including unity 
feedback. 

Then, 




GO 



and the feedback path is H(s) where 



H(s) 



N K s^ + N^K s . 
2 a It 



(n-5) 



(IT-6) 



Therefore, G 2 (s)H(s) is found from Eqs. (II-5) and (II-6) as: 

N2,KaS(s+^^) (II.7) 

^ +60 ** 

The location of the system zero at 8 - — will fluctuate during 

any given cycle of the system and the gain, given by bON^K^, will also 
fluctuate. The difficulties encountered in further root-locus analysis 
are obvious. 

However, there is one useful piece of information which may be glean- 
ed from the root-locus plot. One root of the characteristic equation will 
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nlv/r^ys He betv;cen 5 ^ s ^ 62,016. Therefore, even the nonlinear 
system will alv/ays anproximate a second order system since the other two 
roots w i 1 1 b e d orn i. n a n t . 



The system of F\e, II-l was simulated on a Donner Scientific Corpora- 
tion Analog Computer, Model 3100. The analog computer simulation is 
shown in Fig. IJ-3, The 0 signal is picked off before the final inte- 
gration steo. The simulation for 0 was accomplished after 0 was 
derived as shovm in Eqs. (IT-8). 

9(5) /lO z 

Y(s) “ 

S0(s)/io = 2 V(s) - 2 9(s)/lo 
9(s)/l0 = 2Y(s) -ze(^)/io 



By using the above method for simulating © , differentiating circuits, 
which tend to saturate with steep wavefronts, were avoided. Much of the 

work involved the use of initial conditions. The initial condition for 

« 

0 was set first. Then the initial condition of the integrating ampli- 
fier, whose outmic Is the signal marked "Y" in Fig. II-3, was adjusted 
so that the sumi'w'np amplifier for 0 had the desired initial value of 
Q as its output. Since this system does not represent any parti- 
cular physical system, the terms 0 , and its derivatives, inputs, and 

error signals will be used with units in volts. 

The experimental procedure to be followed in this thesis will be 
to look at the characteristic equation using Mitrovic's method in con- 
junction with the l >2 ~ plane. After predictions about system per- 

formance are made, the v)iedictions will be tested using the analog computer 
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FIGURE H-3 

AMAL06 COMPUTER SIMULATION, (^Rescstarxe in megohms^ ccnpa.Li\anLe. u\ mitvofava^Ls^ pot $et -for sa+uvatifin voltage) 














1 



'h' 






slnulr> tion . Tho result:?; of: computor runs v;ill bo rsnalyzod and compared 
to the pi'GrUctions which v/ero made. Finally, an attempt will be made 
to expand any successful methods of predicting performance for applica- 
tion to other systems X'/ith two nonlinear elements. 
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CHAPTt: ( III 



ITROVIC'S METHOD APPLIED TO THi-: CONTROL SYSTEM 
WITH TWO CAIN-VARIA3LE NONLINEARITIHS 
For the feedback cositrol system described In Chanter II, the charac- 
teristic equation may be nut in the form: 



A,S^ * B,S^ ♦ ■» A„ = 0. 



(iii-i) 



The rppropi'iete Hitrovic^s equation pair from Appendix B is: 

" A,w’0(s>] 

" -wT [■ ♦ A} <0 J'S'>] . 

From ]:!qs. (III-2) it is evident that the parametric equations for 3^ 
and are not functions of the nonlinearities but have constant co- 

efficients which are ~ 1 and A = 60. Substituting these constants 
and the annropriato v'( ^ ) functions from Appendix B in Eqs. (III-2) 
results in (’ll -3). 






(III-3) 



To study this system for stability only, a mapping of the equtals 
xero line of the s-planc into the B^ - plane is desired. iJIicn ^ 

equal s zoro , 



and 






n - ^ 

'2 * 



(III-4c) 



(III-4b) 



can be eliminated from Eqs. (111-4). Thus, the parametric equation 
for stability becomes 
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#• 



1 









"l®2 ■ 



<III-5) 



When values of and v;hich satisfy Eq. (III~5) arc plotted on 

the ~ piano, a graphical solution for the stability curve is 

obtained. This is shown in Fig. III-l. Also shown in Fig. III-l is the 
curve obtained for ^ equal to 0.5. 

The values of and from the characteristic equation are 



given by Eqs. (I I 1-6), 



B, = NiKt 



= 3 ffeON^K^ 



(ni-6) 



A linear M point may be defined as a point on the plane v;hen 

the magnitude of the nonlinear gains, N.j and N 2 in Eqs. (III-6), are 

unity, i.e., when tlie system is linear. 

llith and both equal to one, the linear M-point, from Eqs. 

(Ill 6), is at equal to b2 and equal to 63. There are minimum 

values for and xvdrich occur when N.K or N^K become zero. 

1 2 1 t 2 a 

These minimum values are f^lvon by Eqs. (Ill -6) as = 2 and 3. 

These minimum regions are lined out on the 3^ - 3^ plane in Fig. III-l. 
Since and can have maximum values of unity, then = 62 and 

“ 63 define an upper limit for M-point movement. These boundaries are 
sho\rci as dashed lines in Fig. III-l. Therefore, movement of the M-point 
is completely bounded in a region where 2 ™ B^ = 62 and 3 = = 63. 

LIMIT CYCj..KS. There are no limit cycles for this system by the reasons 
x/hich follov7. Eqs. (III-6) may be rex^ritten as: 
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1 




(III-7) 



B, ^ z 



GO ( Esati 

161 



6=5 + 60 I £ Sail 

^ l*e*l 



If the M'point raoves into the unstable regicn for a significant period 

of time, the oscillations in the system will tend to grov , This growth 

• • 

resi\lts in larger peak magnitudes of acceleration, 0 , and velocity, 

0 . As the magnitudes of 0 and 0 increase, 3^ and become 

smaller, according to grus, (IH-7), x/nich drives the H-point farther into 
the o.'istc'^ble region and the neak raagnitudes of tliO oscillations increase 
even nnre. C<^nverscly, if the M-point does not move into the unstable 
re^'^ion, the peak magnitudes of any oscillations X70uld decrease and the 
A point would move in the direction of the linear M~point. In theory a 
Limit cycle could exist » Hox^jevar, it would be an unstable limit cycle and 
the slightest deviation from such a limit cycle xvould create either a 
stable or an unstable system. Tlierefore, in the actual system no limit 
cycles can exist a.n^ rone v;erc found. 



M pniMT dOTlOp^ T\ie :i-ooinc does not move in a linear system. Therefore, 
ooeration of the system is defined by the M-point location, tx;o roots are 
specified (s = - 1 “ * )> the third root can be found, 

and system performance can be predicted. The M-point location on the 
3^ - nlane can be adjusted to give any desired value of ^ and Ol)^^ 

by adjusting the acceleration and velocity feedback gains. 

If only the tachonetor feedback clianncl saturates, M-point motion 
is a'^ong ;; hori^.ontal line since can vary but is fixed by the 

v^lue of acceleration feedback gain. Conversely, if only the accelera- 
tion feedback channel saturates, movement of the M-point would be along 
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A vertical path. 

I'/hou both feedback paths saturate, the M-point can move anywhere 
in the '• elane within the boundaries previously defined. The 

exact nature of the It point motion can be predicted only if the exact 

• • • 

0 vs 0 ohase portrait can be predicted. However, a qualitative 
description of "■[ point motion can be made using only a few approximations. 
Consider the transfer function of the linear system: 

9CS) _ 60 

e^cs) S5 -.-63$^ H-62S ’ 

which may be rearranged as 

0(s) (s* - 63 ♦ 6Z S ♦ 60) = 60 8i(S) . (ni-S) 

Ati anproxima tion e[ (III-B) is made by factoring out s -i- 62.016 

from the left side of the equation and discarding this term. Then, Eq. 
becomes : 

©CS-)(S^- .184$ .181)= feOe,(s). (ni-9> 

For zero input to the system, Eq. (ITT -9) becomes 

0(s')(s^ ,‘^84s ■+• eS8^) - 0 . (iii-io) 

Ec. (III-IC) may be chanpod to an. approximate differential equation for 
the s y s t em ^ 7b j c t. j ; 

0 -»* 0 0 = 0 (iii-ii) 

IKthin tbio tolerances of the approximations made, Eq. (IIl-ll) will 
govern the linear system. If the system is driven by initial conditions, 
then, at t = 0+, Eq. (lll-ll) must be satisfied. Also, at t = Of, 0 
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hrs not hnc). tine Co chnnf;c ^ron its initial value. If 0(0) is zero, 
then, 

0 =-0 forO<t<CH-^ (III-12) 

To test the validity of the assumptions which have been made, the linear 

system was checked out on the analog computer. Fig, III-2 is a phase 

• » • 

portrait of the 0vs9 plane. In the brief time when 0 t < 0+, the 
approximation given by Eo, (III-12) applies. This time period is indicat- 
ed by the 24 single trace lines in Fig. III-2, Regardless of the Initial 
value of acceleration, when 0(0) ~ 35 volts, 0(0+) = -35 volts. After 
t - 0-f, each of the 12 rims for positive values of 9(0) have essentially 
the same phase portrait. Similarly, each of the 12 runs for negative 

values of 9(0) have the same phase portrait after t = 0+. (It is 

3 

incerescing to note here Chat, if the s tenn were eliminated from Eq. 
(III-8), no significant difference in the ensuing Eqs. (III-IO, 11, and 
12) would have resulted. Use of this fact is made at a later point in 
this cha;ter,) 

Consider this system operating in an underdamped condition. The 
system is then oscillatory. If the velocity has a sinusoidal waveshape, 
the acceleration has a vsimisoidal waveshape which will lead the velocity 
xvv^veshapG by 90 degrees. When the velocity is at a positive or nega- 
tive maximum, the acceleration v;ill be zero. 

If the criteria developed through the aid of the assumptions in the 
preceding paragraphs are applied to the ^ qualitative 

description of M-point motion can be made. Assume the analysis of M- 
ooinC motion is begun vjhen the velocity passes through a maximum value. 

The acceleration is, therefore, zero and the acceleration feedback channel 
is linear. Figure 3 shows the location of this point as . When 
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Phase portrait of the *0 vs S plane^ linear system. 




5TABIUITY CuRve 




i Figure III-3. Predicted M-point motion on the B- - B_ plane. 



the nr?fi:nitMdc of the r^cce Toro t ion bccones f^reater than the saturation 
voltage the magnitude of decreases according to Eq. (III-7) and the 

M noint moves downward as shovm by the arrov; from in Fig. Ill -3. 

This motion \ culd continue until the velocity begins to change its magni- 
tude, say at When this happens, the velocity feedback amplifier 

saci!ratevS to a lesser degree and begins to increase. The M-point 

is then shown by the arrov; from M 2 * A short time later the 
velocity channel goes out of saturation and the acceleration channel 
rercches a maximum degree of saturation. The M-point has then moved to 

As the velocity channel is driven into saturation in the opposite 
direction and the acceleration channel becomes saturated to a lesser de- 
gree, the M-point moves av;ay from in the direction indicated, Assum 

ing a damned system, the M-noint could then oass through M, and arrive 

' ' 4 

at . If, at the acceleration channel remains unsacurated for an 

anpreciable amount of time, the M-point would move in the indicated direc 
tion as the magnitude of the velocity signal decreases. At the ac- 

celeration channel saturates again and the entire process is repeated 
through M- to Mq . At Mq the velocity channel remains unsaturated 
as the acceleration magnitude decreases to This process continues 

until the linear M-point is reached, after which neither channel satur- 
atevS, The above discussion has been for a theoretical M-point motion. 
Before taking up actual M point movement, the theoretical aspects of 
stability will be discussed. 

STARTJJTY 

Theoretically, if the system begins operation in the stable region 
of the plane, there is some damping ratio, , w^hich will 
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cr.nsG the \yeten to have successively str.aller peak signals at the maxi- 
mum of each complete cycle or half-cycle. Therefore a stable system 
would remain stable and M -point motion would follow a trajectory similar 
to that shown in Fig. III-3. Conversely, if the system begins operation 
in the unstable region, then a growth in the peak magnitudes of each 
sxjcceeding cycle occurs and the system can only remain unstablee However 
from Fig. 111-4, a set of initial conditions can be imposed on the system 
so as to place the M-point in the unstable region, say at An un- 

stable system may be cyclic as \ 7 ell as unstable. For a cyclic system 
driven with initial conditions, Sand © will pass through zero. As 
either one masses through zero the nonlinearity imposed on the system by 
saturation is remeved. On the plane, if the value of either 

parameter is toe same as the value at the linear M-point, the instantan- 
eous M-point in question must lie in the stable region. (M^) . Addition- 

ally, at any instant when is greater than 20 ( 0 < 10 volts) or 

• * 

is greater than 30 ( 0<6.66 volts) the instantaneous M-point must 
lie inside the stable region ’ 

Therefore, f^r the non-linear system, the stability curve Is not an 
in' iolate boundary! it muse be possible to start in the unstable region 
and pass into the stable region. Two questions then arise. 

1. Can a system be started in the stable region and become 

unstable? 

2. Can a system be initially unstable and go stable? 

Equation (III-IO) was obtained by assuming that the system at hand 

was essentially a second order system. As previously stated, the same 

3 

result would have been obtained if the s factor had been dropped from 

3 

En. (1II-8). Equation (III-13) results when the s factor Is disregard 
ed but the non-linear coefficients are retained. 
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(III-13) 



OCs) [(3 +(^0N2 ,KJs^ -►(Z-^60NiK^')S + go] = 0, 

Thus, the non-linear differential equation is approximated by* 



(5 + &0Nj,KjS * (Z +foONi.Kt)6 +60e £0. (m-14) 

To move the M-noint to a desired location as a starting point ^ initial 

« • • 

conditions are imrosed on 0 and 0 while 0(0) ~ 0. Thus^ Eqo 

pusy he re^^ritten with the values of the nonlinear gains inserted and for 

the time period 0 < t < Or as. 



(’•T5t)S 



-15) 



for E - 3 volts in each channel, 

sat 

Consider Initial conditions of 0 ( 0 ) - -h 30 volts and 0 ( 0 ) - + 

10 volts. This condition places the M-point at in Fig« If 

F{q, (III-15) is rearranged as Eq , (III-16) and solved for 0 when 0= +30 
♦ * • • 
volts then. 0“ -IrtO volts and the coefficient for 0^ will be at 

4,29 shotm as 

§ = - 0(o) Aqn 0(o)J ciii-16) 

* * 

In order for 0 (Of) to arrive at -140 volts^ it must pass through zero. 
Therefore, M-point motion would be through to and the system 

would be driven immediately into the unstable region. 

Consider the same point but with 0 (0) ~ -30 volts. Then the 

solution tn Eq . /I1I-16) is 0 ~ +20 volts and the M-polnt of t ~ 0+ 
would bo driven to without passing through 

If the system is started at by having 0 ( 0 ) +10 volts and 

• • 

© (0) - +60 volts, where would the M-point be at t = 0+? AgatRj, Eq. 
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fill lb) IS solve'! for 0 i'iie solution is 0 =^ -126 .it'- .''oi , tl^.ore- 

£ore+ 4.43. The M-poimt would move to passing through 

Tor 0((^) ^ -10 volts snri 0 0 :i) ^ -foO volts, then 0 (0-f) w’ouM be ±b . 1 
volts ( ^ 30). This seioence drives the M-point to M . without pass- 

ing through M, . . 

Although the above discussion, is based on approximations seme idea 

of inliTial fP pcint movement has been obtained. When initial conditions 

« «• 

are inpO'sed on d and 0 so that they both have the same sign^ the system 
apparently would be hit harder than if B and ^ were of opposite sign. 
Thus, the question of stability would hinge^ lo part^ on the question of 

the sign of the initial conditions imposed o 

*0 

When £nu is solved for 0 under varying initial, conditions 

0 » 0 

for 0 and 0 * labie results. In lieu of all the apprcxi'mat ions 

whicl'j have been made^ Table III-l can riot be regarded as an exact predic- 
tion of system stability. It does indicate the type of initial condition 
settings which would tend to make the system stable^ marginally stable^ 
or comf/ietely unstable. At this point It is concluded that e'xperimental 
tests must be 'conducted to determine stability limits for this systerrio 

To investigate stability^ the system was driven with various initial 
conditions.. The starting noirkt on cha plane is defined by the 

A #4. A A 

magnitudes of the 0 and Q initial conditions. The 9 vs 0 phase 

nlane was selected as the appropriate method of recording data because 
of the direct correlation beween this phase plane and the 
plane. 

After conducting several computer runs 5 it became obvious that soiree 

• # * 

regions of the 9 vs 0 plane would produce stable runs and some regions 
w'ould prod'uco an unstable systenu To delineate the stable regior.;,, the 
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TABLE III-l 



©(0) 


8(0) 


%(o) 


BgCo) 

♦ 


' Initially: 
Stable (S) 
Unstable (U) 


8(0*) 


B2(0+) 


Predicted: 
Stable (S) 

Iftis table (U) 
Marginal (S/U) 


+15 


+30 


lU 


9 


S 


-130 


li.38 


s/u 


+15 


+6o 


Hi 


6 


s 


*-130 


Ii.38 


s/u 


+l5 


+90 


Hi 


,5, 


s 


-130 


li.38 


s/u 


+15 


-30 


Hi 


9 


\S , ■■ 


- 10^ 


. 21.00 


S 


+15 


-60 


Hi 


6 • 


s 


-10 


21.00 


S 


•!-+15 


-90 


Hi 


■5 


s 


- 10 


21.00 


S 


-15 


+30 


Hi 


9 


s 


+il0 _ 


21.00 


S 

i. 


-i5 


+6o 


Hi 


6 


s 


* + 10 


21.00 


S 


-15 


+90 


Hi 


5 


s 


+;io 


21.00 


s 


-15 


-30 


Hi, 


9 - 


s 


+130 


.'li.38 


s/u 


-15 


-60 


Hi 


6 


s 


+130 


li.38' 


s/u 


-15 


-90 


Hi 


5 


s 


+130 


li .38 


s/u , • 

\ 


+30 


+30 


9 


9 




-HiO 


li .29 


u . 


+30 


+6o 


9 


6 


u 


H 

1 


li.29 


u . 


+30 


+90. 


9 


5 


u 


-lliO 


li.29 ■ 


U ; 


+30 


-30 


9 


^ 9 


. 'S 


- 20 


12.00 


1 s 


+30 


-60 


9 


6 


V'’ 

■'U 


- 20 


12.00 


s • - 


+30 


-90 


9 


5 


u ^ 


- 20 


12.00 


s 


+30 


+30 


9 




si' 


+ 20 


12,00 


s 


-30 


+6o 


9 


6 


u 


+ 20 


12,00 


s ; 


-30 


+90 


9 


5 


u 


+ 20 


12.00 


S 


-30 


+30 


9 


9 ' 


s 


+HiO 


li .29 


U 


-30 


-6o 


5 


6 


u 


+HiO 


li .29 


u 


-30 


-90 


9 


5 


u 


•+HiO 


li .29 


u 
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Q viy Q nlanc was .searched until a stable run lesulited from one 

initial condition but^ an unstable run resulted from an initial condition 
disniaced a small dist»ance from the original starting point. Figures III- 
ua through IIT-bg are n sampling of these type runs. VJ^hen the iocus of all 
stable starting points is drar;n, the 0 vs © phase plane is divided 
into a stable region and an unstable regiono This is shown In Fig. III- 
7, For ail these runs, 0 (0) equals zero* 
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i.II 7 tuo Idea thr.t stability depends on the signs 

of else itpti.n.l cond-itlnns as well as the magnitudes. When the locus of 
stab I o ia;ns is nlocied on th.c plane, as In Figure III-8, two 

^•rinaipai div^idinw lines result. The first, slio\^i as a dashed line, is 

* A • 

roi initial conditions of 0 and 0 of opposite signs. The second 
^ivi ting line results from & and 0 intiial conditions of the same 
s J gn . 

fin. axi-titativc prediction of M-point motion is confirmed, in part, 
i;hev) ihii mns of Fig. are ?>l»atted on the plane, ns in 

i'p 111 p -re run is nor predictable from the qualitative des- 
cription n W r. 'itir.n. 

* • 4 

>'lll t ■ 0 -J:i © trajectory can be predicted, stability can not 

« • • 

bo rredic lew To assist in predicting the ^ vs ^ trajectory, isocline 
d'leory seens to be the only available method. The differential equation 
for this control system is rewritten as; 

§ + B^6 + 6^0 <o0 0 =-0. (III-17) 



ii . 

dt “ ^ 



<III-18) 







dt 




(ITI-19) 



t: n suhstitiaing os. (IilI-IB and 19) in Eq . (ITI-17) and rearranging 



gi ves ; 




dt 



B^0 -B^9 -600. 



(III-20) 
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T» ^ 
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T 






I'k 0 



0 



ohasc nlnne i s cU? mod as N, then 



N i il. 

d.6 

When Che lefc side of En. ( 111 - 20 ) is divided b> 

« • 

side by 0 , !•>. (Ill 22) results. 



( 111 - 21 ) 

iXBfdt and the right 



de/dt 

de/it 



Eliminating dt 


irorn cite 


by die slope,, N, 


produces 



_ ^ B g, 0 ^ Bt G (qO 0 

““ # s 

0 



( 111 - 22 ) 



left side of EQc (1X1-22) and replacing dS/d^ 
the desired equation for the isoclines in the 



0 vs 0 rhase plane When the non-linear factors, and are sub- 

stituted In the isocline equation, it becomes: 




- 60 © 



( 111 - 23 ) 



Louptioa Iii-23 rney be rewritten as: 



N — -30 -ZQ - Co 0 -180st<tr\ 6 - X8o (m- 24 ) 



• « « 

Ihc only vr?kr.nvn^ v?hen Eq. (III-24) is applied to the © vs 0 phase 
plane for given values of H, is © . If some value of 0 is assumedj the 
isoclines can be dravm. Figures Ill-iOa through 10c show these isoclines 
for 0 enual to zero, 10* and 20, respectively. 

From iig. HI 7, it appears that, for a stable system, the maximum 

• S 0 

vaUie of 0 when © is zero is approximately 20 volts. By starting with 

« 

© at -20 volts in Fig. Ill -10b* and tracing a trajectory away from that 
noira in a reverse oirection, a path results sxich that, for initial condi- 
tions anywhere on the vath and for 0 remaining at a value of +10 volts, 
the trajectory viould return to the starting point. For a positive velocity 
rnd a positive displacement, the displacement can only increase. Thus, 
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« t • 

t'\c -rr ■ t O'- jrc tory in tlic unper half of the 0 vs 0 plane v;ould 

t^'v/rrd .Mjo 'U-'i incd fho i.r.nclines in Fiy,. III-lOc. For S equal 
t:<‘ */ero, • r>econa preFiccoc! liraicctory to arrive at that point is dravm 
on Fig, 111 K'b. Using these two predicted trajectories, a region of 
stable runs is predicted or. Fig. Ill-lOb. 

Once 9 v';0 nhase plane was searched for stable and 

rrv.u.c’oFrr regions voer 0^ > ^^'^s ten volts. Figures Ill-lla through llg 
rrc a sc'iin» :-Tr\ ^0 '-tins .ii id) were Tna.de in this search. When the locus 

(O' st'lilc starrinp •'oj.nts is drawn, the re.gion of stable runs is enclosed. 

Lg' re Tj I '? sb. n- s n -•'rtion of this region vdrich corresponds closely 
v; i 1 1 i t ' 1 o r e , i ' n v:l \ i c h as 'U e d i. c t e cl . 

^'ccording t:o the results polotted in Fig. III-12, vdicn 0(0) is 
d'U) volts no rositive vadues of 0 are permitted if the system is Co 
re.’ain ^tablc . '[iierefore. it is concluded that the system cannot remain 
sLrh]o. for a sten treater thaT> 10 volts by the reasoning which 

foViov/s. Since t'v"* .'sto • f s i>ot critically damped, for a step input of 
10 \ there • il -=vre overshoot. As 0 reaches 10 volts a positive 

• elcclt.) t : : t d-vc overshoot. Therefore, the system will go un- 

st^'^ble when © flO vo^tc- and the velocity is positive. 

Figure III ]3a ia a vr;rse portrait of the 0 vs 0 plane showing the 
tr^^lectory for a steo input ^\f 10.1 volts. For this input, the system 

rc'’>ri*ncJ stabh, 'p ir L: 1 3b siu>ws n recorder trace for the indicated 

« 

V- ,ci '*bl es verrns t ’ . "i' rc ill Id s the abase portrait of the 0 vs 0 

p’oa a nlrne .b. i t 1 .: 'vv.it vas iucreasod to 10.2 volts. For the 10,2 

volt in nit, tin unstable. No limit cycle operation could 

be hy vT' 'f . luuit between 10.1 and 10,2 volts. When the input 

is tSuG in sr; '• stops, less than ten volts, and the system has 
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l] '• to rtv ivCc-ndy state between each su'cessv^^e step input the 
- iM V>'\v^ to ony nagtiitudc dosired. 
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rive chosen to place the linear M-ooint on the stability cejrtc, an un- 
stable Unit cycle results, An example of this type of run is shown in 
FJy. Ill IF for the cor responding point from Fig, 111-17, 

Ircra/^sinp the iJecdljack gains K tind K should widen the stable 

^ a t 

• t • 

royiion ox the 0 vs 0 phase plane since the M-point would spend a long- 
ei: neriod of time in the stable region of the - B 2 plane. Converse- 

ly, reducing the feedback gain ; should reduce the stable reglom These 
tv’o nredic t: i oas were found co be true. Figure III-19 shows an unstable 
run r-o^‘uits T 7 bcn is reduced to 0,13 whereas, with ~ 1«0, the 

st'^rtiu'^^ ’oint xras In the stable region. Figure ITl-20 shows two stable 

nms rdien V ^nd K are increased to 10. Both of these runs were un- 
a t 

stable n’itii unity feedback gains. 



It is concluded that an accurate prediction of system performance 
cannot be made using Mitrovic’s method in conjunction x>/ith algebraic and 
graph?Cr?1 methods vrhen tx/o coefficients contain non-linear factors. From 
the 1 n\^est i " OTIS undoi taken by this author, it appears chat isocline 
tb.ooiy by <tsc>lf v;ou 1 d provide the most accurate method of prediction)., 
Ik'wever, tho use of x isocline t-'eory v/ould require a large number of 
isocline nlots, with accuracy being proportional to Che number emplcyedo 
An analysis x^ 7 ould proceed in the follox^ing manner, given the initial 
conditions at the starting noint. Using the following relationship^ 

and the values of ©, ^ short trajectory 9 compute At, 
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select a ncu i^>c. lir.c ri -t i.cccrdir 2 g lo the nev va.i;- .A & » : rr; ^ 

second shore traject.:>r> on the net'' .Iso.:, liar, plot pro,' o. a.- t . 

This process would hi' repeated until the ttajectory eirh-r . *':> t" - 

co'Ties unstable. The isocline theory developed here Vvceld I- 

cable to systems uUli i bar at t e r is t Ic equations cf higher orhec tnan 
third order equrt C" U>r this system 



It is also concluded that the b-est means available <-.’,^.1 ai- any 

system v^ith tv;o non-1 Ineai i ties is by analog or digital c .on . r str :y 

Unfortunately^ computer methods can answer the question o: vn-ij hap:- ■ 

but not the more important quest ion of why. 

The dualism ot the equals zervO line on ibe T, - h :ia.;c is Ihe 

i - 

single most disturhing arc’ i’".<-yqvl i ah i - phLromenon f :i In A 
If one such duaiisti: lire exists it stems entirely pcssiM' that ocn^rs 
might also he lonrJ ..1ms laetor al'-v^ne would preciudt- th-, ot 

Mitrovic'’s neinofi to pr--.dict system per fonriance Tvrhen twv:« '/k n- 1 ro^.d i t ies 

are involved. 
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CHAFFER I¥ 



M1TR0¥IC"S METHOD APPLIED TO THE STABLE CONTRO^u SYSTEM 
In Chapter Illg the control system was studied to ascertain stability 
characteristics of the control system when the gain was high encugh to 
cause the uncompensated system to be inherently unstable « By reducing the 
gain, the system can be made inherently stable. From Eq, which 

is repeated here as Eq, (IV-l)^ the parametric equation for stability isz 

- K, (I¥~l) 



which may be rew'ritten as 

(2 -f KN K )0 + KN,K ) K, (If- 2) 

It' 2 a 

When K and K are zero, the maximum value which K may have fer a stable 
at 

system is six, which is found from Eq. (I¥-2). The character! site equa- 
tion for the system becomes 



+ B + B..S + 6 = 0 
1 ^ 



The solution of Eq, (l¥-3) with the appropriate B^-E^ Mitrcvic equation 
pair from Appendix B gives the desired parametric equations for the Bj^- 
plane: 



n 






6 ( 2 ^ ) 



(W'-4) 



6 + 






n 
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The graphical solutions to EqSo (IV-4) are shocm in Figo IV-1* Sin^e 
cannot be less than two and cannot be less than threej the corres- 
ponding areas are ruled out in Fig. !¥-!• ThuSj the M-point must alvays ba 
in the stable region and the control system is inherently stable. No dual- 
ism of the zeta equals zero line could be found by varying either initial 
conditions or the size of the step input. 

It is desirable to predict an average M-point motion on the 
If such a prediction were possible , system performance could be accurately 
described, Ti^o approaches to the prediction problem were attempted. The 
first approach was a linear approximation to a step input response. The 
second approach was to analyze the results of many runs to ascertain If 
some definite *'pattern" of performance could be determined, 

LINEAR APPROXIMATION METHOD, 

The procedure in this method was to choose a step input for the system 

and a starting M-point. It s assumed throughout that the system could be 

regarded as a second order system and^ therefore, that the values of zeta 

and natural frequency, would govern system performance. 

The step input size was chosen as ten volts. Initial values of - 

2.57 and =0.6 were chosen to place the starting M-point at = 17,4 

and = 6.6. At this point, the value of zeta is 0,4 and O - 3,8 rad/sec, 
z n 

If 0 were sinusoidal then 9 would lag 9 by 90 degrees and 9 would lag 

9 by 90 degrees. Therefore, if 9 goes from its initial value of zero to 

0 

the maximum value, M , in time, t, then 0 would arrive at its maximum 

pt 

• # 

value in t/2 and 0 would have a maximum in time t/4. 

Defining 9 as, 

^ max ’ 
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~ 9(0 ) 



0 - 9 

max nvax 






and 9 as 
max 



9 - 9 

max max 



(I¥»6) 



t/2 

provides a basis for predicting M-point motion on the plane. 

For zeta equal to 0,4 the maximum overshoot would be lo25<, Therefore , 

for a ten volt inputs, ® is 12,5 volts. For 6J equal to 3,8 rad/sec^ t 

max n § 

« 

equals 0,414, Substituting these values in Eq. (IV-5) gives a of 

30,2 volts. Therefore. . equals 2.2?1, The values of 9 is prudict- 

ed from Eq. (IV-6) as 146 volts j which corresponds to a of 3,206, 

Thus, the M-point would be bounded by the straight dashed lines shown In 

Fig, IV-2. A revision to the first approximation is now made by choosing 

an average trajectory along the zeta equals 0,2 line and guessing at an 

average value of 4) of 1,43 rad/sec, from the corner value of 0 , Under 

n n 

these conditions, M is 1.5 and t equals 1,10 seconds. Using this value 

pt 

of M ,0 is 15 volts and, from Eq, (I¥-5) , 9 is 13,6 volts and^ there- 
pt max max 

A A 

fore B, . is 2,52* The revised value of 0 is computed from Eq, CIV-6) 

1 min max 

as 24,3 volts, which corresponds to a . of 4,23, A second revision 

2 man 

is now made based on the curved dash line of Fig. IV using a value of 

c 

zeta equals 0,3 and an 4) of 1,4, From this data, M is 1,359 9 is 

a mi>: 

12,0 volts and 0 is 21,3 volts. Corresponding to “ 3,50 and 

max ^ ^ I 

. - 3,85, Sirce this second revision shows about the same M-point 

2 min ^ 

motion (shown by the curved dashed line in Fig, l¥-2) no further revisions 

are made. 

To predict M-point motion after the first peaks are reached ^ it is 
assumed that in one complete cycle positiouj velocity, and acceleration 
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are all damped by a factor oi € 



The follovir^.a d' 



point motion ovoir any one cycle is made assuming that and E. non- 

linear throughout the entire cycle andj thereforej, ha^'ie val^i-cs gi'eo 07 

B. ^ 2 + 18 



and 

B, - 3 + 18 , 

' 1¥T 

It is also ass'unied that velocity and acceleratiofii have an a-^^erage 

tude over one cycle equal to the average magnitude of a sine vave o^^er cae 

cycle ,, i .6 , 



and 



Therefore.^ 



0 



ave 



,636 0 



peak 



,636 0 , . 
ave peak 



11 ave 



^ 2+18 



and 



12 ave 



2 + 

>636 10,, j 



A B “ B ^ - B 

1 ave 12 ave 11 ave 



,636 



1 G 1 


2 S 


1 1 


1 2« 




1 1 1 



•1 



then 



AB, = 18 

1 



.636 |9J 

where 6 is ths damping factor. 



1 - £ 
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Similarly 



, B ^ 28 o 4 ^ 
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The S faccorg appearing in Eqs. (IV-7 and lV-8) are in 

Table TV-1. Using the values of S from the table and Eqs. (1/-7 ana 
lV-8), a predicted average M-point trajectory is plotted in iEig. IV-3. 



^ TABLE IV -1 



0.0 


0.00 


0.1 


0,105 


0.2 


0.222 


0.3 


0.352 


0.4 


0,493 


0.5 


0 640 


0.6 


0,823 


0.7 


1.011 


0.8 


1.222 


0.9 


1,461 


1,0 


1,719 



The results of the computations for average M-point trajoctoty are sho'vn xn 
Table IV-2. 

The revised predicted trajectory for a ten volt step Input was based 
on peak magnitudes of velocity and acceleration. If a single point average 
value is used based on the average value of a sine viave. over a complete 





2 


I oj(yolts') 


TABLE IV-2 
1 9 j volts) 


s 






3.50 


3,50 


12.0 




36.0 


.105 


.25 


,08 


3.75 


3.58 


10.3 




31.0 


.110 


.30 


.10 


4.05 


3.68 


8,8 




26.5 


.140 


,45 


.15 


4.50 


3,83 


7.2 




21,7 


,165 


.65 


o & 


5.15 


4,05 


5,7 




17.1 


.210 


1,05 


,35 


6.20 


4,40 


4,3 




12.9 


.285 


1.88 


, 63 


8,08 


5.03 


3,0 




8.9 


.423 


4,00 


1.35 


12,08 


6.38 


1.8 




5.3 


1.719 


27,10 


9.20 


*17,^0 

^ is 

** 6 ^ is : 


•>^*6 5 60 

linear for 
linear for 


B, - 17.4 , 
1 

= 6,6 . 
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cycle chon 6 q.63 Q and 0 = 0.G3Q* . For the t r ' .iit step 

avG 

# a » 

input* 0 “ 7.)5 volts and 0 - 13.4 volts correspondini^ t: a 4 

^ avG ave l a . e 

= 4,39 and a B. == 4«13. This point is shown as IL la Fig lV'-3o 

2 ave ^ 1 ave ^ 

According to the predictions developed here^, if the first M is as 

i av6 

shown in Fig, I¥-3 then about three complete oscillations should take place 
before the average M-point arrives at the linear M-point. 

Figure l¥-5 is an M-point trajectory for a 10 volt seep for 

• •• 

which the 9 vs 0 phase portrait is shown in Fig, I¥-4o Is only one 

complete oscillation in the actual system. Figures IV-6 through l'v-i2 are 
phase portraits and M-point trajectories for step inputs of 20 ^ 30 ^ 40 5 and 
50 volts. 

The average M-point trajectory does follow the general trend cf the pre- 
dicted average trajectory, Howeverj the predicted trajectory allows for 
nearly twice the number of cycles as actually occur. 
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CON^n:SiONS. 



Intre axe primary sources for error in the predicted M-pcjr 
trajcv^cory. The first source is in predicting the initial of exCc- 

ity and acceleration and is due to the delay in position tespcnse after 
the step input is applied. This results in erroneous values oi tfie fltsc 
veU-: ity and acceleration peak magnitudes* 

iJie second source is in assuming a sinusoidal response vh^n tt?:e actual 
response is a damped sinusoid ♦ The damped sinusoid has a lover average 
magnitude over a cycle than apure sinusoid and, thus, the piedicted rumber 
of t^x^piete os,.iilations are greater than the actual nurriber 
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(.AS T . 



'ihc 



conpu ter 
t Tc>n s Te r 



:c ’T^ -T. ^ Y 
U VTT. ';! Til A 
!> I n c ! '. d ■' r r r . i.* 

ninulntlon of the 
function for this 



•TTii TiTii nriini/: iSTiijn -/ 
Sl'YOU C1!/^I:^TL. 
i:iiG system to ho snnly^.cd 
system are shown in Fig. 
system is 



and the nnalog 
A-1. The closed 



1 oop 



- 10 N 

QjlS) s’ * 3 - 2S * ION 

where M represents the instantaneous variable gain of the saturation 
nonl inerrlty . I>y choosing the last tv?o coefficients of the characteris- 
tic G^^uatlon as the variable coefficients, the characteristic equation 
be cones 



33^ ^ B^S ^ = 0. 



(A-2) 



From Appendix B, Hitrovic’s equations for 3^ and are: 

B„ = - 

and 

= 300^0^1:5) * <o^03(^) . 

I^'ubsci tuting the values of the 0 functions for ^ = 0 from Anpendix B in 
Fes. (A~3) and (A-4) gives the parametric equations of the stability 
c u r VO a s 

= 36 jJ rt- 5 ) 

and 

. (A-6) 
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\7Lt^ vrhios of 



• T. 



c'md 



r ( • > •^•'rlod rr> 6J^ ?n v.'^riod, Tliur, tlic ranbilitv 

cur\-*- 1 slwvn. ’y conpr r . (A-'l) ^n. t ') (A 2), ? c> r.ooii lhni‘ 

?j , ™ 2 = ION where K has a "-laxinun value of one. Thus 

an ?! ^o;int V)cus is specified on the - 13^ piano. 

'P'lc intorsoction r>f the fi noint: locus with die stability curve defines 
a sjndc '‘oint xf'icro all die e'^uati.ons for and ^ are satisfied. 

At iAv'f' '^inc, 



= 3 coA = 10 N (A- 7 ) 

and 

= 6d^ - Z . fA-^:) 

Tiu.^ so’uti ni o,r : . / 'Ives " 1.414 rad/soc, Suhsti tutinp, this 

vahio O)^ a. (A-7) yiodhs N “ a . ^'inco the variable fyain of 

the '^aturativi non ’ i noa r i : v is dc'inod as the rat.o of the output to the 
ll'p - r n. Lu'b> • tiio error siynal can be predicted fron Ka. (A~9) 



P ^ T T - r , , ^ .; V ^ ^ i . 



£ 




VO 1 1 51 



(A^9) 



The linear system ’litb a vain of 10 is unstable, Tiius, if any 
disturbance ' cciirs in tlie system, the eutnut signal v;ill begin to in- 
crease "in i-aynitude. Th uv^Pore, tiie error signal also increases in 
nagnitiuhm ‘.an the r>rror channel saturates the gain of the system is 

decreased. <- tiie rovMO'^ gain i5", decreased, tlie system becomes more 

•Atrhie and A. rror sjvuia} begins to decrease. A decrease in the error 
signal increases the system gain and the above process is repeated. 
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"vcM'itiia Mv r Avnrnic i ] ihrin n is rcnchcfl in \;hich the systi i osci . ":it‘. s 

alniv'; t’.o :! -point locus at the intersection of the H-point lor us /-rut the 
stability curve* iiue value o], the nonlinear pain^ 13, an'l the 
vof the error sirnal^ S , arc therefore tine averaged values » 

Fig. A~3 is 'I nhase portrait of this system obtained with an analog 
cenputer si''ul/>tion ni the system. The radian frequency of lo395 rad/ 
sec. agrees ^;ichin one percent of the predicted value of 6)^ = 1.414 rad/ 
sec. As previously stated, the magnitude of the error signal^ £ . when 
the system is operating in a limit cycle can he predicted as 8.33 volts, 
f^^hen the error channel saturates, the system will be driven with a con- 
stant input oi five volts Therefore, the system should be driven at the 
sene rare in either a positive or a negative direction and the error sign- 
al channel should snend the same amount of time in positive saturation as 
it does in pepative saturation. The predicted magnitude of the error 
signal should then be an average value which is independent of the sign 
of the error signal. On Fig. A-3 a vertical line was drawn at £ - -8v4 
volts. Then, for the limit cycle, the area marked ’’A'* can be compared to 
the area marked by counting the number of squares contained in eacli 

area. The two areas are equal and, therefore, the average magnitude of 
the error signal is .,4 volts which agrees closely with the predicted 
VP lue . 

Fig. 3-4 shows the graphical solution to this nonlinear problem on 
the 8^ plane and the *■ plane. The solution for the Ih-, ~ 

plane follows the same pattern as that for the plants 

Solving this problem on the plane is, however^ somewhat unique 

and is explained in the following paragraphs. 
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.'r(v iV'C'!):: ’ , L‘H :’Ltrov.ic oqxiPMono foi" this sysicn 

rpnropr) /)ic c oc i c \(Mit :'ron the characteristic equation ir-.s. ■ ‘-d: are 

and 

^7- “ jj” ^ 

IHicn the ) functions are substituted in Eqs. (A~101) ani 'A-ii '^hey 

be cone ^ 

Bo = w'] 

and 

®2. " (z^) [^ “ 6)^ (1 ’ i 

The suhstitution c' ‘ iu Eqs, (A-i2) and (A-13) tesulrs In singulari- 
ties in both L'^uat'ont^ Hut ^ for 0, s = j 6 l)^ . Sutsti '! -ting this 

'ujantity for s in fa Af-2) yields 

'jK + zjQ„ + =0. (4-iv. 

Requiring the real and imaginary parts of Eq. (A-14) to g: to 
independent ly produces the following equation paiCc 

jcoJ7-io^) = 0 

Bo - = 0 . fa-*6 

o 

The solution to >q. (;\-I5). after first dividing out j is 40^^^ 

Thus, 60^“ rad/sec and from Eq ,- CA-lb)^ - 2E^ Si:. • f. - 
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3 j.rom the characteristic equation, then = 6 and N = 0,6, Zii. 

- B 2 curves are then plotted and a graphical solution is obtained 

v;hich is consistent with the solutions from the - B, and Bvv “ f- 

0 1 0 i 

planes . 



CASE II. SATOILATION IN THE VELOCITY FEEDBACK CHANNEL. 

The block diagram of the system to be analyzed and the analog 
computer simulation of the system are shown in Fig. A-5o The closed 
loop transfer function for this system is 



6, (a _ 30 

0;(a s'* . T 5" I- 14 + (8 *-iaON}s * 30 






where N represents the variable gain of the nonlineari tyo As in Case 
I, the last tv;o coefficients of the characteristic equation are chosen 
as the variable coefficients. The characteristic equation is 



3"^ -*■ 7s^ + + B^s + 6^ = 0 



(A-18) 



where 



and 



B. = 30, 



B^= 8-^lZON. 



From Appendix Mitrovic’s equations for and B are! 

u 1 



Bo-- 






(A-19) 

(a-20) 



«'A-2i) 



and 



, A k,, 
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FC'.'dba.- k control syjtein with ^ aVu-.-ar i-'n xr. the t v<homot«-r ice.lbark chann 



Substitutiu," Che */■( ) funccions for '^=0 from AppendiA 5 ir. Eqs 

mid (A-22) provider. Clio following equations for the stability cjives, 

Bq ~ 14" Ct)^ ■“ . (A-i3> 

The stability curve and M-point locus are plotted on the P- - h 

U I 

on Fig, A''6„ 

The graphical solution to this problem reveals both a stable ard an 
unstable limit cycle. The gain of the linear system (N - l.G) is high 
enough so that the 1 incar system is unstable. Any distuirfeardce that occurs 
in the system producers an oscillatory state, UTien the velocity feedback 
signal Sc?turaces the feedback channel, the value of the nonlinear gain is 
reduced. The M-point moves to the left toward the stability cutve. If 
the average M-point moves to the left of (or inside) the stability curve 
the system is damned. As the oscillations decrease in magrditade the nor - 
linear gain increases rnd the M-point returns to the stability curve. 
Eventually a dynrnic equilibrium is reached and a stable limit cycle re- 
sults, At the intersections of the M-point locus and the stability 
curve Eq, (A 19) and Eq. (A-23) are equivalent and^ therefore^ 

- 14 -*■ 30 - 0. (A-25) 



The solutions to Eq. (A-25) are: 



and 



60^ = 1.G3 v'cid/sec 



6)^- 3. 37 rod/sec. 



/A -26) 



(fA-27;* 
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At each Interscc t ion of the H-point locus Eqs j A-». 

are c^;ui valent ami 

120 N = - 8. A-.8 

Pot the stable Imit cycle using the value of f ro'in Eq. ■( A-x 7 •' , eAe 
solution ot fA-cH) is N equals 0.596o For the unstable iiait ■. y l- 

. (A~26) applies and the solution to Eqo (A-28) is N equals ''’bP8 'u, 

Ihe unstable limit cycle results when the magnituie •.^1 t':-‘ 
feedback signal is large enough to produce saturation su;t: that ^ is 
equal to or less than 0.0874^ This condition will exist loir a-'-- v rragn.- 
tude of 0 equal to or greater than 5^73 rad/seCc, a value 
the following equation: 

N = -,t I < 0 . 0874 . 

6|e I 

Foi ,^ny 0 whose magnitude satisfies (A-29)^ the systtm is - 
stable and the oscillations begin to grow. This Increase in “h" magroilL.de 
of 0 further reduces the value of N and the oscilliati''^'ns : 'cri. rw- t r 
greater :nagnitudes. Fig. A-7 is a phase portrait ot this system shcA*. na, 
operation at both limit cycles. 

It ii> emphasized here tliat the operation of the system at rkiths’r 

limit cycle Is nor in a small region about the intersection >.^1 the 

« 

stability cut'/o and the M-point locuSo As the value of © passes through 
zero in either direction^, the value of N Is one^ since, tne vel^wity i-el- 
back channel is rot sataiatedo ThuSg the instarttainieous valcss ci 
cause the system to operate through a long portion of 

The unst'd.>1e limit cycle^, for example^ would produce values ji F - 

JL 

tv/cen eight end 128, 
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Lh?>r m sAit.i.An(;w w the acceleration feedback channe:,. 

Th'? block of the system to be analyzed ard the 

com[<nter siiiiul at i on of the system are shown in Fig« A-8o Th- cI.'S- 
loop transfer function of tlic system was found from whi-ih the cr ara 
istic equation is" 

s'* + 7s^ + (I4 ^10N)s^ + 8 5 20 = 0 



where , 



and 



= 8 

= 14 +10N. 



[k-ji ^ 



■:'A- 



The Mitrovic equation pair are selected fccm Apr'-;di,n 

With the correct coefficients from Eq* (A‘-30) and the 0( ^ lun-.-'t iacs 
for ^ Oj the parametric equations for the stability cerve axes 

B;i = 7coJ A-'o) 




^4 j) 



At the intersection of the H~point locus and the stability carv-e £qs. 
(A-31) and (A-33) are equivalent« Therefore^ ~ f 

this value of 63^ ^ equals 18o64j from Eq, (A-34) ^ and equals 

0,4645 from Eq, (A-32) , 



Ine graphical solution to this problem is sho'wr. in Fig. If 

the magnitude of the acceleration signal is large enough sc-ch that A' 

Is less than 0,46a_^ then the. system is unstable and osol 1 lat i o‘:jS la- 
crease in magnitude, This increases the peak magnitude of th-:;' accviera- 
tlon sipnrd , red’irces the magnitude of and causes the syst--m t- - 

ome mere unstable. Therefore the system will never attaiLW a stable 



For small signal magnitudes of acceleration the system op- rates nt-Jr 










i' 

'j 





ontrol syscem with saturation in a»' ; - ie r ar loo ft"t:dba^k Lhann 



I 





■tSfe 



the point wliere i: ic unity* In this region the system has sen^^ dair^.- 

in,g ratio, , greater than zero. As the oscillations are decreased in 
magnitude^ the value of N approaches unity. In this mode of eperatlen< 
the oscillations v>?ill eventually die out and the system will operate vn 
a linear^^ stable state. Fig. A-10 is a phase portrait of this system 
which illustrates the stable mode of operation^ the unstable limit cy:le 
and the unstable mode of operation. 
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table li. Equatioas for the 0j^( ^ ) fuactions aFp,ar 

'4 ft) = 0 
( 0 = -1 

4 = 2 p 

%(\) = l- 4 f 
if ft) = "4f 

0sOi>^ -lfe'5'' 

6 ^ - 32 ^^ + 32 ^^ 

0,ft)= 1 - 24^^ +80^'' - 6+5^ 

« 
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